Far-field emission profiles from L3 photonic crystal cavity modes 
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We experimentally characterize the spatial far-field emission profiles for the two lowest confined 
modes of a photonic crystal cavity of the L3 type, finding a good agreement with FDTD simula- 
tions. We then link the far-field profiles to relevant features of the cavity mode near-fields, using 
a simple Fabry-Perot resonator model. The effect of disorder on far-field cavity profiles is clarified 
through comparison between experiments and simulations. These results can be useful for emission 
engineering from active centers embedded in the cavity. 
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I. INTRODUCTION 

Photonic crystals (PhC) offer unprecedented control 
over electromagnetic field confinement in all three spa- 
tial directions In particular, two-dimensional PhC 
nanocavities in a planar waveguide have already found 
applications in different fields such as nanolasers, non- 
linear optics and quantum information processing 0, Ej[ . 
Similar to any electromagnetic resonator, PhC nanocav- 
ity modes are essentially characterized by two figures of 
merit: the cavity quality factor, Q, and the effective con- 
finement volume of each mode, V mo< \e (3- The quality 
factor is proportional to the photon lifetime in the cavity 
which depends on the cavity losses to the external world. 
The mode volume is a quantitative measure of the spa- 
tial confinement of the electromagnetic mode. In most 
applications, it is crucial to maximize the Q/V mo d e ratio. 
For example, the Purcell factor, which measures the en- 
hancement of the spontaneous emission rates for atoms 
resonant with a cavity is directly proportional to this fig- 
ure of merit. In PhC nanocavities the mode is strongly 
confined to a very small volume, on the order of (A/n) 3 , 
where A is the mode wavelength. In a planar membrane 
nanocavity, in-plane confinement is provided by spatial 
localization of a structural defect in a perfectly periodic 
PhC with a photonic band-gap, while out-of-plane con- 
finement is given by total internal reflection between the 
slab and the air cladding (assuming a suspended mem- 
brane as a planar waveguide). Very high quality factors, 
in the range 10 4 -10 6 have been demonstrated in the 
literature. In particular, the L3-type cavity, consisting of 
three missing holes in a triangular lattice, was the first 
PhC cavity to show quality factors larger than 10 4 @. 

The spectral mode structure for L3 cavities has been 



thoroughly investigated by Chalcraft and coworkers [8| 
who compared the calculated resonant energies, quality 
factors and emission polarizations for the lowest-order 
modes with experimental data. Most experiments cou- 
pling single quantum dot emitters to a nanocavity exploit 
the fundamental (i.e. lowest-energy) cavity mode fgL fTol[ . 
However, higher order modes can still be important for, 
e.g., efficient pumping in nanocavity lasers [11], selective 
excitation of quantum dots embedded within the cavity 
fl2l [l3j , or mutually coupling quantum dots in different 



spatial positions [141 ] . Several groups have studied the 
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near-field emission profiles of photonic crystal nanocav- 
ities [Til . fToj ] . even with polarization- resolving imaging 

In this paper we report an experimental and theoretical 
investigation of the spatial far-field profile of the out-of- 
plane emission for the two lowest-order modes of L3-type 
PhC nanocavities. We believe the characterization of the 
out-of-plane far-field emission for PhCs is important for 
two main reasons. First, for single-photon source appli- 
cations the emitted radiation needs to be efficiently col- 
lected into a fiber, and simultaneous optimization of far- 
field emission for multiple nanocavity resonances could 
be useful. In addition, in the case of cavity-QED exper- 
iments in the "one-dimensional atom" approximation, a 
perfect mode-matching is needed to get a large enough in- 
terference between the input light field and the field radi- 
ated by the atom [l8l - |20j . Recently, quite some work has 
been done to get a beam-like vertical emission from PhC 
nanocavities [2l|-[24j for the fundamental mode. Here we 
extend previous work by experimentally analyzing the 
far-field emission properties of both the fundamental and 
the second-order mode, finding good agreement with nu- 
merical simulations. We introduce a simple model, based 
on a one-dimensional Fabry-Perot resonator, to estimate 
the essential far-field characteristics of a given near-field 
mode profile and link them to the relevant device param- 
eters. We believe that such a model can be useful for 
fast parameter optimization, while full-scale numerical 
simulations can provide an accurate but time-consuming 
description of the electromagnetic field in the structure. 
Finally, we will discuss the effect of fabrication impcr- 
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fections on the far-field cavity emission profiles. As we 
will show, measurements of far-field profiles are relatively 
easy to perform and they can provide insightful informa- 
tion about the parameters and the quality of the cavities 
under examination. 

The paper is organized as follows: in Section II, we 
present experimental measurements of the far-field pro- 
files and a comparison with theoretical far-fields ex- 
tracted from finite-difference-time-domain (FDTD) sim- 
ulations; in Section III, we introduce a simple model, 
based on a Fabry-Pcrot resonator, which is sufficient to 
give indications of what the actual far-field profile looks 
like for a given near-field and to link far-field properties 
to actual device parameters. 



II. THEORETICAL MODELING AND 
EXPERIMENTAL DATA 

Our sample consists of a 180 nm GaAs membrane 
grown by molecular beam epitaxy on top of a 0.92 /mi 
Alo.7Gao.3As sacrificial layer on a GaAs substrate. An 
Ino.4Gao.6As quantum dot layer is grown at the center 
of the GaAs membrane by depositing 10 periods of 0.55 
A-thick InAs and 1.2 A-thick Ino.13Gao.s7As. The L3 
PhC cavities were fabricated on the sample using stan- 
dard electron beam lithography and reactive ion etching 
techniques [25|, [26| . The lattice constant of the triangu- 
lar hole lattice is a = 240 nm. The L3 cavity design was 
properly modified for Q optimization (see modified holes 
in Fig. P H3,[23. 

The sample was placed in a He-flow cryostat at about 
5 K and illuminated above the GaAs bandgap with a few 
mW laser beam (wavelength 780 nm) on a few /j,m 2 spot. 
The photolumincsccncc from the quantum dot layer em- 
bedded in the membrane was collected in the direction 
normal to the membrane using a microscope objective 
with numerical aperture NA = 0.8 and spectrally ana- 
lyzed with a spectrometer (resolution 5.5 GHz/pixel). An 
example of the spectral emission is shown in Fig. [T] for 
a device with R = 54 nm. According to theoretical pre- 
dictions based on a guided-mode expansion method [29| , 
the cavity supports two confined modes, with resonances 
respectively at A^ = 987 nm (theoretical Q-factor 
Qth ~ 180000) and A^ = 957 nm (Q th ~ 15000). Ex- 
perimentally, we measured the first-order mode with a 
Lorentzian profile centered around Ai = 982.5 nm with 
a full-width at half-maximum (FWHM) of 0.195 ± 0.002 
nm, from which we extract an experimental quality factor 
of Q ~ 5000. On the other hand, the second-order mode 
has a less perfect Lorentzian lineshape centered around 
A 2 = 956.4 nm with FWHM 0.63 ± 0.03 nm, and an 
experimental Q-factor Q ~ 1500. Experimental quality 
factors are lower then the predicted ones due to scatter- 
ing from fabrication imperfections (30j and possible ab- 
sorption from sub-bandgap trap levels and surface states 

Experimentally, the emitted radiation collected from 
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FIG. 1: Top: sketch of the L3 photonic crystal cavity used in 
the experiments. Given a design radius R for the holes, the 
two holes closest to and in-line with the cavity have reduced 
design radii R' = 0.75R [2^] . The small holes are shifted 
from their lattice positions out from the center of the cavity 
by a fixed quantity (55 nm). Finally, the third holes (i.e. two 
holes away from the small holes) are shifted from their lat- 
tice positions out from the center of the cavity by 45 nm [28j], 
to further increase the theoretical Q. Bottom: experimen- 
tally measured photoluminescence spectra for the two cavity 
modes. 

both modes results in a strongly linearly polarized sig- 
nal. However, as it is shown in the 3D FDTD simula- 
tions of Fig. [2 each near-field mode profile has x- and 
y-components of the electric field of comparable inten- 
sity. The reason for the detection of linearly polarized 
light can be found by calculating the far-field projections 
of such polarization-resolved near-field profiles. 
The far-field profile can be obtained from the near-field 
using the procedure introduced by Vuckovic and cowork- 
ers |32j |. According to the surface equivalence theorem, 
all the information about the far-field profile can be ob- 
tained from equivalent electric and magnetic currents, 
J s = n x H and M s = n x E, which depend on the 
in-plane near field components: 

N x = -FT 2 (H y ) N y = FT 2 (H x ) . . 

L x = FT 2 (E y ) L y = -FT 2 (E X ) (L > 

where FT 2 denoted the two-dimensional Fourier trans- 
form. These equivalent currents are used to calculate 
the retarded vector potential of the electromagnetic field, 
which in the far-field can be related to Fourier transforms 
of the near-fields. The radiation intensity per unit solid 
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angle can be calculated as: 
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where 77 = \J fio/eo is the impedance of free-space and A 
is the mode wavelength. The radiation vectors in spher- 
ical coordinates can be expressed from their cartesian 
components as: 

Ng = (N% cos (p + Ny sin (p) cos 6 , , 

N v = —N x sin ip + N y cos ip , ^ ' 

and similarly for L. The far-field profiles calculated from 
the near-fields in the left panels of Fig. [2] are shown in 
the same figure, on the right. In these plots, the color 
scale is normalized to the totally emitted power in the 
upper half-space of the PhC cavity (the same normaliza- 
tion factor is used for E x and E y ). Most of the emission 
from the x-polarization of both modes is predicted at 
very large angles, and therefore is inefficiently collected 
by commonly employed microscope objectives. This re- 
sults in the strong linear polarization observed in the 
photoluminescence spectra. 

To perform a direct measurement of the far-field emit- 
ted intensity at each resonant mode frequency, the fil- 
tered photoluminescence at the back focal plane of the 
microscope objective is imaged. Given a characteristic 
size of the near-field emission, p ~ 500 nm, at a wave- 
length A ~ 1 fim the Fresnel number is F = p 2 /(LA) ~ 
0.01, well in the far-field regime (L ~ 2 mm). The far- 
field was imaged on an intensified CCD camera by a lens 
with focal length 40 cm in a 2f — 2/ configuration. To 
make sure that we were looking at the microscope ob- 
jective back focal plane, we adjusted the lens to see the 
sharp image of the objective edge on the CCD. This sharp 
edge was used to calibrate the numerical aperture scale 
of the far-field images, assuming that the sharp edges 
correspond to the NA = 0.8 of the objective employed. 
An interference filter, with a bandwidth 1 nm, was used 
to spectrally select the mode of interest. Images were 
collected after integrating for 30 s and the background 
noise was removed by subtracting an image taken with a 
slightly tilted interference filter. 

The experimental far-field spatial emission profiles 
for the first-order and second-order modes are shown 
in the two larger plots on the right side of Fig. [31 
together with the far-field projections obtained from 
FDTD simulations (smaller insets on the left) for a 
direct comparison. The first-order mode exhibits a 
centrally illuminated area extending to about NA ~ 0.5, 
with a ring-like structure inside (NA ~ 0.2), matching 
the low-NA portion of the simulated far-field. The 
simulated far-field suggests that most of the light from 
the first-order mode is emitted in two high-NA lobes, 
which are not collected at all by our set-up. The far-field 
profile for the second-order mode consists of two lobes, 
whose center is at a minimal NA ~ 0.3. 
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FIG. 2: FDTD simulations of polarization resolved near-field 
patterns (left panels) and far-field emission profiles (right pan- 
els) for the first and second modes of an L3 cavity. Each mode 
has two linearly-polarized components. The white concentric 
rings correspond to a grid with numerical-aperture step of 0.2. 
For far-field plots, the color scale is normalized to the totally 
emitted power in the upper half-space of the PhC cavity (the 
same normalization factor is used for E x and E v ). For the x- 
polarization, most of the far-field emission is at large angles, 
making collection less efficient than for y-polarization and re- 
sulting in a strongly polarized photoluminescence signal. 



Finer details also appear inside the two lobes, in the 
form of spots separated by Ak/k ~ 0.1 — 0.2. Such 
structures correspond, in the near-field, to light sources 
which arc separated from the cavity about 5 — 10 times 
the characteristic size of the cavity mode. The simu- 
lated near-field profiles shown in Fig. [2] suggest that the 
optical field decays very fast out of the cavity region, 
implying that such features might be due to light that 
escapes from the cavity due to fabrication imperfections 
[3(il ]. The fine details are reproducible for different mea- 
surements performed on the same device. Finally, we 
show in Fig. 2] the far-field profiles of the second-order 
modes for different devices on the same wafer. Each plot 
shows the characteristic two-lobes profile, as expected for 
this mode, with reproducible finer details that appear to 
be device-dependent. 
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FIG. 3: Experimental far-field emission profiles for the two 
L3 cavity modes (larger panels) and the corresponding pre- 
dictions by FDTD simulations (smaller panels on the top left 
of the experimental plots) , resulting from the sum of x- and y- 
components in Fig. 2. The white concentric rings correspond 
to a grid with numerical-aperture step of 0.2. The color scale 
bars show the detected photon counts per second. The counts 
for the two profiles cannot be compared since they depend 
not only on the device parameters, but also on the density of 
quantum dots emitting in the wavelength range of the specific 
cavity mode. 



To investigate the cause of fine structure within the 
two-lobe far-field pattern of the second order mode, 
FDTD simulations were done in which disorder was in- 
troduced. It is assumed that, due to fabrication imper- 
fections, the dominant type of disorder is in the hole radii 
of the PhC lattice. Therefore, the hole radii R were 
varied randomly according to the distribution function 
P(R) oc exp[— (R — R) 2 /2a 2 )}. Some simulation data are 
given in FigureEJ showing disorder introduced to the far- 
field profiles as a result of increasing the disorder parame- 




FIG. 4: Experimental far-field profiles of the second-order 
mode emission from four different devices on the same wafer. 
Finer details are reproducible in different measurements per- 
formed on the same device, but differ from device to de- 
vice. The white concentric rings correspond to a grid with 
numerical- aperture step of 0.2. 
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FIG. 5: Far-field intensity profiles for varying degrees of dis- 
order, normalized by the total power radiated into the upper 
half space. The images are cut-off at a numerical aperture 
value of 0.8, indicated by the white line, to allow for better 
comparison with experimental data. Data are shown for dis- 
order parameters of a = 1 nm (left column), 2 nm (middle 
column), and 5 nm (right column). 



tcr a. These results indicate that far-field measurements 
could be used as an indicator of disorder in the lattice 
structure. 
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III. A SIMPLE FABRY-PEROT MODEL 

In general, 3D FDTD simulations can provide accurate 
modeling of near- field and far-field properties of PhC cav- 
ity modes. However, they give little physical insight on 
how the detected features of such far-field profiles can 
be related to specific device parameters. In this Section, 
we will show that the experimental data can be repro- 
duced by a simple model, elaborated from the proposal 
of Sauvan and coworkers (33j . 

For a line of N > 3 missing holes, the PhC nanocav- 
ity can be described quite accurately by a Fabry-Perot 
resonator, in which the fundamental Bloch mode of a 
single-linc-defcct PhC waveguide is trapped between two 
PhC mirrors of modal reflectivity r(X). The properties 
of such a cavity are shown to depend only on three pa- 
rameters, namely the group index n g of the Bloch mode, 
the reflection coefficient r(A) of the mirrors and the effec- 
tive cavity length L. The Bloch mode can be calculated 
as the eigenstate of the PhC waveguide in the Fourier 
basis, and its modal reflectivit y c an be obtained with 
the method described in Ref. [3J]. Fabry-Perot mod- 
els have been shown to be a useful tool to probe cav- 
ity resonances and the group index of photonic crystal 
waveguides [H, [36| , and to describe acousto-mechanical 
cavity tuning effects [37|. We consider the E y near-field 
profiles shown in Fig. [2] Taking the intensity distribu- 
tion along the y = axis, the modes show a sinusoidal 
intensity distribution with nodes and antinodes in the 
cavity region, decaying exponentially outside the cavity 
region. From simulations, the intensity decay length for 
the first-order mode is Si ~ 220 ± 10 mn, while it is 
5-2 ~ 155 ± 8 mn for the second-order mode. This is 
very similar to the intensity distribution in a cavity be- 
tween two Distributed Bragg Reflectors (DBR). For the 
E y near-field profiles most of the radiation is emitted 
along the y = axis, with little structure and weaker 
intensity outside. Given the relevance of such polariza- 
tion in determining the far-field emission properties, to 
a first-order approximation it makes sense to consider 
a one-dimensional model, which takes into account only 
the structure along the y = axis. Such a model has the 
advantage of being extremely simple, although able to 
give significant hints on the main far-field profile prop- 
erties. A two-dimensional model based on an effective 
index approximation should be used (with no significant 
difficulties) to study also the E x profiles. 

Let us consider a one-dimensional Fabry-Perot res- 
onator, with the field intensity profile varying sinu- 
soidally in the cavity region and exponentially decay- 
ing in the mirror regions, with characteristic penetra- 
tion depth S. The resonant frequency £l m can be calcu- 
lated by imposing the total phase acquired by the Bloch 
mode traveling back and forth to be a multiple of ir. 

The far-field profile can be calculated using the proce- Eq. H Let's consider a separable electric field distribu- 
dure outlined in the previous Section (Eq. Since tion E y {x,y) = a(x)fi(y). The component L x isjep- 
the modes are TE-like modes, only E x , E y and H z are arable as well: L. x — a(k x )(3(k y ), where a and j3 are 
non-negligible and only the L x component is relevant in the one-dimensional Fourier transforms of a(x) and /3(y), 




FIG. 6: Far-field results for the one-dimensional Fabry-Perot 
model (first and second-order modes). For each mode, on the 
left side the expected near-field profile considering a Fabry- 
Perot cavity with DBR mirrors with varying penetration 
depth 5. On the right side, the corresponding far-field profiles 
calculated from Eq. [4] For the first-order mode, in the region 
\k x \/ko < 0.5, there is a single peak centered at k x /ko = for 
small penetration depth. This peak flattens and broadens on 
increasing penetration depth. For large enough penetration 
depth, the central peak splits into two smaller peaks. The 
second-order mode, in the region | fcn, | / fco < 0.5, consists of 
two peaks that split further and further out for increasing 
penetration depth 5. 
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FIG. 7: One-dimensional far-field profiles as a function of 
the intensity penetration depths calculated using the simple 
Fabry-Perot model for the first (left-side plot) and the second 
(right-side plot) order modes. The white dashed lines corre- 
spond to the experimental position of the center of the lobes 
(Fig. O and the penetration depth from FDTD simulations. 

making: 

K(9,ip) oc |a(fc o sin0cos</2)| 2 |/3(fcosin#sin!y9)| 2 x r.\ 
x (cos 2 ip cos 2 9 + sin 2 tp) 

In a simplified onc-dimcnsional model, j3(y) is narrow 
in real-space, so its Fourier transform is wide and can 
be neglected {\(3(k y )\ 2 ~ 1). If we look at the far-field 
distribution along the x-axis we select ip = 0, so that the 
resulting one-dimensional far-field profile is: 

K(6) ~ |a(fc sin<9)| 2 cos 2 6> (5) 

or, in terms of transverse wavevectors: K(k x ) ~ 
\a{k x )\\kl-kl). 

An analytical solution can be found and is reported 
in Appendix A. However, the resulting formula is too 
complicated to give intuitive insights; therefore we will 
just discuss the numerical results (from Eq. |12[) in this 
context. Fig. [5] shows the expected far-field profiles for 



the first and second-order modes for different values of 
the penetration depth, 8. The first-order mode exhibits 
a central peak (centered at k = 0) and two outer lobes 
as predicted by the FDTD simulations. We see that, for 
increasing penetration depth, the two main outer lobes 
become narrower in k-space and more outward, while the 
central peak broadens and flattens. A similar behavior 
can be observed in the second-order mode. Here, the 
far-field profile is given by two central peaks and two 
outer lobes, and again, for increasing penetration depth 
8 the two outer lobes move outwards and become more 
localized. 

The Fabry-Perot resonator model was shown to give 
results for the quality factor Q i n ag reement with more 
sophisticated FDTD simulations [331 ] : 

Q = Y^ n 9 ( L + 2L p) ( 6 ) 

where the group delay r = experienced by the light 
upon mirror reflection enters as a characteristic length 
L p = — Xq/ (Aim g )(d<l)r / d\)\ . In general, for PhC struc- 
tures (for example, heterostructure mirrors) it has been 
shown that L p can be unrelated to the characteristic 
damping length of the energy distribution inside the mir- 
rors 8 [381 ] . However, in the configuration under investi- 
gation, the classical relation for a DBR, L p < 8, is valid, 
with the equality being strictly fulfilled only in the limit- 
ing case of quarter- wave mirrors with low refractive index 
contrast, An/n (L p = 8 = nA/(2An)). 

In Fig. [JJ the emission in the region \k x /ko\ < 0.5 is 
shown as a function of the penetration depth. In the case 
of the fundamental mode, there is a single peak centered 
around k = for small penetration depth, with a quite 
flat profile and broadening for \k x /ko\ < 0.25. When the 
penetration depth increases beyond 8 ~ 150 nm, the far- 
field emission splits into two lobes, which separate more 
and more on increasing 8. At 8 ~ 250 nm, the peaks are 
centered around k ~ ±0.25fco, corresponding to an open- 
ing angle of about 15 degrees. These findings explain the 
ring-like structure we observe in the experimental far- 
field: the vertical dashed white line in the figure shows 
the experimental central NA for the peaks (from Fig. [3]) , 
which corresponds to a penetration depth of about 220 
nm, fully compatible with the penetration depth from 
FDTD simulations, 8i. For the second-order mode, in 
the region (fc^/fcol < 0.5 there are always two emission 
lobes, even at small penetration depth. On increasing 
8, the two lobes move further and further apart. From 
Fig. [3J the lobes are centered around NA ~ 0.35, which 
in our simulations is consistent with a penetration depth 
of about 150 nm, as predicted by FDTD simulations. 
Therefore, far-field profiles can provide useful informa- 
tion about the penetration depth 8 of PhC cavity modes, 
which in this case offer a bound on the effective cavity 
length, L p , and therefore on the Q-factor. Extensions of 
the present analysis to treat coupled cavity modes [39l . l40T | 
can also be considered. 
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IV. CONCLUSIONS 

In conclusion, we have presented an extensive charac- 
terization of far-field emission profiles from L3-type pho- 
tonic crystal nanocavities, introducing a simple imaging 
technique as an efficient tool to give a two-dimensional 
mapping of the emitted intensity. The measurements 
have been directly compared to theoretically modeled 
far-field projections from the 3D FDTD near-field cavity 
modes profiles, and we believe these results to be useful 
to PhC cavity designs for specific purposes. The effect 
of disorder on the far-field profiles was investigated via 
numerical simulations. 

Finally, we have introduced a simple Fabry-Perot model 
that is able to capture the essential features of far-field 
properties for suitably designed near-field profiles. As 
a particular application of this model, we can envision, 
for example, the simultaneous optimization of in- and 
out-coupling for two different modes supported by the 
same PhC cavity, which is still an open problem that 
might benefit from simplified models like the one pre- 
sented here. 
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Appendix A: Analytical calculations 

The resonance frequencies for the modes can be found 
by setting the condition that the total phase acquired by 
the Bloch mode traveling back and forth is a multiple of 
7r (2k m L = rmr), which results in: 



mir 
~2L 



(7) 



Let us start with a perfectly confined mode, with no pen- 
etration into the mirrors (5 = 0). In this case the field is 
given by: 



a s=0 {x) = n 



where Ii{x) is the rectangular function, Ii(x) = 1 for 
|x| < 1/2 and zero elsewhere. The phase <j) m is set by 
the boundary conditions: <p m = for m even (cosine-like 
solutions) and 4> m = tt/2 for m odd (sine-like solutions). 
From Eq. |U 
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(9) 



[*™ ( S in0-1)] + 

+e- , '* m Sine [*f- (sin(9 + 1)] | 2 cos 2 
The two modulo-squared Sine functions in Eq. [9] give 
two main peaks centered at 9 = ±7r/2, which correspond 
to the higher-NA peaks in the far-field in Fig. 6. The first 
relative maximum of Sinc 2 (a;) is at x ~ 1.43037T, which 
for Eq. [9] corresponds to sin# ~ ±(1 — 2. 861 /to). The 
FWHM of such a peak for Sinc 2 (x) is 0.522-7T, which cor- 
responds to Sk x ~ 1.044/to (Sk x ~ 0.35 for the first-order 
mode, corresponding to to = 3, and 5k x ~ 0.26 for the 
second-order mode, corresponding to m = 4). Therefore 
for the first-order mode, the first relative maxima of the 
Sine functions superpose, giving just one central peak. 
For the secod-order mode, on the other hand, the first 
relative maxima are well separated. 

Including the penetration depth S into the model, a 
simple near-field profile can be taken as a superposition of 
II and two exponentially-decaying wings, as follows: 

a{x) = [n (f ) +H(x-%) e -(*-V2)/* + 

+H (-x - |) e -(-*- L / 2 )/ 5 ] e <(*m*+*m) + c. c. 

(10) 

where H{x) is the Heaviside function (H(x) = 1 for x > 
0, H(x) = for x < 0). The part around x = (II (f )) 
gives the same Fourier-transform as in Eq. [9] (which we 
label &i(sm0))), while the left and right-side exponential 
decay regions give the following: 
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e im(7t/2)(sine-l) 
S/L— im7r(sin — 1) 

-im(7r/2)( 8 me-l) 



6/ L + 7777 7T (sill — 1) 
(11) 

All the quantities depend on the ratio S/L between the 
penetration depth and the cavity length. The far-field 
profile can be calculated to be: 



Ki(0) 



i=i 



cos 2 9 



(12) 



(8) 
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